Scaling properties of step bunches induced by sublimation and related mechanisms: A unified perspective

This work provides a ground for a quantitative interpretation of experiments on step bunching during sublimation of crystals with a pronounced Ehrlich-Schwoebel (ES) barrier in the regime of weak desorption. A strong step bunching instability takes place when the kinetic length is larger than the average distance between the steps on the vicinal surface. In the opposite limit the instability is weak and step bunching can occur only when the magnitude of step-step repulsion is small. The central result are power law relations of the between the width, the height, and the minimum interstep distance of a bunch. These relations are obtained from a continuum evolution equation for the surface profile, which is derived from the discrete step dynamical equations for. The analysis of the continuum equation reveals the existence of two types of stationary bunch profiles with different scaling properties. Through a mathematical equivalence on the level of the discrete step equations as well as on the continuum level, our results carry over to the problems of step bunching induced by growth with a strong inverse ES effect, and by electromigration in the attachment/detachment limited regime. Thus our work provides support for the existence of universality classes of step bunching instabilities [A. Pimpinelli et al., Phys. Rev. Lett. 88, 206103 (2002)], but some aspects of the universality scenario need to be revised.Comment: 21 pages, 8 figure

Equilibrium of anchored interfaces with quenched disordered growth

The roughening behavior of a one-dimensional interface fluctuating under quenched disorder growth is examined while keeping an anchored boundary. The latter introduces detailed balance conditions which allows for a thorough analysis of equilibrium aspects at both macroscopic and microscopic scales. It is found that the interface roughens linearly with the substrate size only in the vicinity of special disorder realizations. Otherwise, it remains stiff and tilted.Comment: 6 pages, 3 postscript figure

Short-time scaling behavior of growing interfaces

The short-time evolution of a growing interface is studied within the framework of the dynamic renormalization group approach for the Kadar-Parisi-Zhang (KPZ) equation and for an idealized continuum model of molecular beam epitaxy (MBE). The scaling behavior of response and correlation functions is reminiscent of the ``initial slip'' behavior found in purely dissipative critical relaxation (model A) and critical relaxation with conserved order parameter (model B), respectively. Unlike model A the initial slip exponent for the KPZ equation can be expressed by the dynamical exponent z. In 1+1 dimensions, for which z is known exactly, the analytical theory for the KPZ equation is confirmed by a Monte-Carlo simulation of a simple ballistic deposition model. In 2+1 dimensions z is estimated from the short-time evolution of the correlation function.Comment: 27 pages LaTeX with epsf style, 4 figures in eps format, submitted to Phys. Rev.

Drift causes anomalous exponents in growth processes

The effect of a drift term in the presence of fixed boundaries is studied for the one-dimensional Edwards-Wilkinson equation, to reveal a general mechanism that causes a change of exponents for a very broad class of growth processes. This mechanism represents a relevant perturbation and therefore is important for the interpretation of experimental and numerical results. In effect, the mechanism leads to the roughness exponent assuming the same value as the growth exponent. In the case of the Edwards-Wilkinson equation this implies exponents deviating from those expected by dimensional analysis.Comment: 4 pages, 1 figure, REVTeX; accepted for publication in PRL; added note and reference

Driven Lattice Gases with Quenched Disorder: Exact Results and Different Macroscopic Regimes

We study the effect of quenched spatial disorder on the steady states of driven systems of interacting particles. Two sorts of models are studied: disordered drop-push processes and their generalizations, and the disordered asymmetric simple exclusion process. We write down the exact steady-state measure, and consequently a number of physical quantities explicitly, for the drop-push dynamics in any dimensions for arbitrary disorder. We find that three qualitatively different regimes of behaviour are possible in 1-$d$ disordered driven systems. In the Vanishing-Current regime, the steady-state current approaches zero in the thermodynamic limit. A system with a non-zero current can either be in the Homogeneous regime, chracterized by a single macroscopic density, or the Segregated-Density regime, with macroscopic regions of different densities. We comment on certain important constraints to be taken care of in any field theory of disordered systems.Comment: RevTex, 17pages, 18 figures included using psfig.st

Lonely adatoms in space

There is a close relation between the problems of second layer nucleation in epitaxial crystal growth and chemical surface reactions, such as hydrogen recombination, on interstellar dust grains. In both cases standard rate equation analysis has been found to fail because the process takes place in a confined geometry. Using scaling arguments developed in the context of second layer nucleation, I present a simple derivation of the hydrogen recombination rate for small and large grains. I clarify the reasons for the failure of rate equations for small grains, and point out a logarithmic correction to the reaction rate when the reaction is limited by the desorption of hydrogen atoms (the second order reaction regime)

Adaptation dynamics of the quasispecies model

We study the adaptation dynamics of an initially maladapted population evolving via the elementary processes of mutation and selection. The evolution occurs on rugged fitness landscapes which are defined on the multi-dimensional genotypic space and have many local peaks separated by low fitness valleys. We mainly focus on the Eigen's model that describes the deterministic dynamics of an infinite number of self-replicating molecules. In the stationary state, for small mutation rates such a population forms a {\it quasispecies} which consists of the fittest genotype and its closely related mutants. The quasispecies dynamics on rugged fitness landscape follow a punctuated (or step-like) pattern in which a population jumps from a low fitness peak to a higher one, stays there for a considerable time before shifting the peak again and eventually reaches the global maximum of the fitness landscape. We calculate exactly several properties of this dynamical process within a simplified version of the quasispecies model.Comment: Proceedings of Statphys conference at IIT Guwahati, to be published in Praman

Dynamic Scaling in a 2+1 Dimensional Limited Mobility Model of Epitaxial Growth

We study statistical scale invariance and dynamic scaling in a simple solid-on-solid 2+1 - dimensional limited mobility discrete model of nonequilibrium surface growth, which we believe should describe the low temperature kinetic roughening properties of molecular beam epitaxy. The model exhibits long-lived ``transient'' anomalous and multiaffine dynamic scaling properties similar to that found in the corresponding 1+1 - dimensional problem. Using large-scale simulations we obtain the relevant scaling exponents, and compare with continuum theories.Comment: 5 pages, 4 ps figures included, RevTe

Variational Formulation for the KPZ and Related Kinetic Equations

We present a variational formulation for the Kardar-Parisi-Zhang (KPZ) equation that leads to a thermodynamic-like potential for the KPZ as well as for other related kinetic equations. For the KPZ case, with the knowledge of such a potential we prove some global shift invariance properties previously conjectured by other authors. We also show a few results about the form of the stationary probability distribution function for arbitrary dimensions. The procedure used for KPZ was extended in order to derive more general forms of such a functional leading to other nonlinear kinetic equations, as well as cases with density dependent surface tension.Comment: RevTex, 8pgs, double colum

Persistence exponents for fluctuating interfaces

Numerical and analytic results for the exponent \theta describing the decay of the first return probability of an interface to its initial height are obtained for a large class of linear Langevin equations. The models are parametrized by the dynamic roughness exponent \beta, with 0 < \beta < 1; for \beta = 1/2 the time evolution is Markovian. Using simulations of solid-on-solid models, of the discretized continuum equations as well as of the associated zero-dimensional stationary Gaussian process, we address two problems: The return of an initially flat interface, and the return to an initial state with fully developed steady state roughness. The two problems are shown to be governed by different exponents. For the steady state case we point out the equivalence to fractional Brownian motion, which has a return exponent \theta_S = 1 - \beta. The exponent \theta_0 for the flat initial condition appears to be nontrivial. We prove that \theta_0 \to \infty for \beta \to 0, \theta_0 \geq \theta_S for \beta 1/2, and calculate \theta_{0,S} perturbatively to first order in an expansion around the Markovian case \beta = 1/2. Using the exact result \theta_S = 1 - \beta, accurate upper and lower bounds on \theta_0 can be derived which show, in particular, that \theta_0 \geq (1 - \beta)^2/\beta for small \beta.Comment: 12 pages, REVTEX, 6 Postscript figures, needs multicol.sty and epsf.st